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The pulsation equations for spherically symmetric black hole and 
soliton solutions are brought into a standard form. The formu- 
lae apply to a large class of field theoretical matter models and 
can easily be worked out for specific examples. The close rela- 
tion to the energy principle in terms of the second variation of 
the Schwarzschild mass is also established. The use of the gen- 
eral expressions is illustrated for the Einstein- Yang-Mills and the 
Einstein-Skyrme system. 
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1 Introduction 



In recent years the stability properties of numerous new black hole and soliton 
configurations, arising in selfgravitating nonlinear field theories, have been 
studied. Early investigations of this kind were concerned with the instabilities 
of such objects in SU(2) Einstein- Yang-Mills theory. It was shown in Jl|, 
that both the particle-like solutions of Bartnik and McKinnon 0], as well as 
the "colored" non-Abelian black holes [Q are very fragile structures. Some 
black holes with hair which were found in other matter models turned out 
to be (at least linearly) stable. This is, for instance, the case for the Skyrme 
model || D, [7|] or when a Higgs triplet is added to the Yang-Mills fields 
[§, P, |10|. If, however, a Higgs doublet is coupled, black holes and solitons 



turn out to be unstable again [|11 



Most of these and other investigations are restricted to spherically sym- 
metric configurations. A closer look at them reveals many common features. 
Hence, one would like to have a more general treatment, encompassing a suf- 
ficiently large class of field theoretical models. For any given one, it should 
then be possible, for example, to write down readily the final pulsation equa- 
tions, without repeating the numerous intermediate computations. 

A save and systematic method for a stability analysis is to linearize the 
coupled field equations around a given equilibrium solution and to study 
the frequency spectrum of the resulting (multicomponent) eigenvalue prob- 
lem. For spherically symmetric situations one expects intuitively, that the 
resulting fluctuation operator is closely related to the second variation of the 
Schwarzschild mass. Similar "energy principles" are important in other fields 
of physics such as, for instance, in plasma physics for ideal magnetohydro- 
dynamic conditions |T3. In classical mechanics the second variation 5 2 H of 
the Hamiltonian provides the well-known Lagrange-Dirichlet stability crite- 
rion, which applies usually also in infinite dimensions. (A counter example in 
elasticity theory has been discovered by Ball and Marsden [0.) One should, 
however, remember that the energy criterion is often of no use. A well-known 
example is provided by the equilibrium positions of the restricted three-body 
problem. Here one has to linearize the Hamiltonian vector field in order to 
investigate linear stability; 5 2 H is indefinite and gives no information on the 
fluctuation spectrum. 

If the energy principle applies, it provides usually the simplest method 
to decide on (linear) stability. In this paper we establish the validity of 
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the energy principle for a large class of matter models coupled to gravity. 
More precisely: Under certain assumptions on the dimensionally reduced 
form of the matter Lagrangian, we show that the fluctuation operator of 
the pulsation equation can be read off from the second variation S 2 M of the 
Schwarzschild mass. Although this is not surprising, we felt that a general 
demonstration of this is needed, in particular since energy considerations 
in gravity are notoriously problematic. More important, in deriving this 
connection we arrived at explicit expressions for the pulsation equations and 
for S 2 M. These can readily be worked out for a given field theoretical model, 
which should considerably simplify future stability studies. 

The article is organized as follows. In section 2 we specify the form of the 
dimensionally reduced Lagrangian of the matter model. The coupled field 
equations for (time-dependent) spherically symmetric configurations are then 
expressed entirely in terms of this effective Lagrangian, taken for a diagonal 
form of the metric on orbit space (see eqs. ( |2.33| ) - ( f2.36| )). This involves 
some nontrivial considerations since, in the time-dependent case, one would 
otherwise loose one of Einstein's field equations. The general set up is used 
in section 3 to bring the pulsation equations into a standard form, involving 
only the matter perturbations. From these the metric perturbations can be 
obtained algebraically, without solving any differential equations. In section 
4 we demonstrate that the fluctuation operator can be read off from the 
second variation of the total mass. The Einstein- Yang-Mills system and the 
Einstein-Skyrme system are taken as illustrations for ready applications of 
our general formulae. The dimensional reduction of the gravitational La- 
grangian and the derivation of the mass variation formula for spherically 
symmetric configurations are deferred to two appendices. 

2 Dimensional Reduction 

2.1 The effective Lagrangian 

Let us consider a spherically symmetric spacetime (M,g). By definition, 
(M, g) is a Lorentz manifold on which SO (3) acts as an isometry group, such 
that all group orbits are metric two-spheres. This guarantees that spacetime 
(M, g) is locally the warped product of a two-dimensional Lorentz manifold 
(M,g) (the orbit space M/SO(3) = M with induced metric g) and the two- 
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sphere S 2 with standard metric dQ 2 : 

M = MxS 2 , g = g + r 2 dtt 2 . (2.1) 

The positive function r: M — > R (the Schwarzschild coordinate) is assumed 
to have no critical points on M. For the volume form fj on M there is a 
corresponding split: fj — r\ A r 2 dQ, where rj and dfi are the volume forms on 
M and S 2 , respectively. In addition, we introduce the intrinsically defined 
function iV = (dr\dr), where ( • | • ) denotes the inner product on M and the 
mass fraction m is defined by iV = 1 — 2m/ r. 

As is shown in appendix A, the dimensional reduction of the Einstein- 
Hilbert action yields for the effective gravitational Lagrangian L G , after sub- 
tracting the standard boundary term, 

L G V = -jz[dm\dr)rj. (2.2) 

Expanding dm with respect to a coordinate basis, dm = m dt + m' dr, and in- 
troducing the metric functions (3 = —(dt \ dr)/N and S = sj—g, the effective 
Lagrangian takes the form 

L G r] = C G dt Adr , C G = (m' - (3m)S . (2.3) 

In terms of N, S and (3 the metric on M becomes 

g = -NS 2 (dt 2 + 2 (3 dt dr) + (— - (3 2 NS 2 ) dr 2 . (2.4) 

The effective Lagrangian Cq contains, as we shall show below, the entire 
dynamical information. Requiring that the matter Lagrangian depends only 
on the metric and its first derivatives, we shall now establish the fact that 
the Einstein equations agree with the Euler-Lagrange equations for the total 
effective Lagrangian. 

Let us first consider the gravitational part. Defining the Euler-Lagrange 
operator E° f for a dynamical variable /, say, according to 




we immediately find from eq. fl2.3| ) 

E° m C G = -S', E° s C G = m', E°£ G = -mS. (2.6) 
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Using well known expressions for the Einstein tensor G^ v in Schwarzschild 
coordinates, we see that 

E° m C G = -LSN- 1 {G\-G\), (2.7) 

E° S £ G = -y^> (2-8) 

E°£ G = ~SG\. (2.9) 

Let us now turn to the effective matter Lagrangian Cm, defined by 

S M = G- 1 J C M dtAdr , (2.10) 

where Sm denotes the matter action. By virtue of the definition of the energy 
momentum tensor T^, the variation of Sm with respect to the dynamical 
variable / = N, S, (3 yields 



SfS M = ljT« b { d -§f5f}^ 



df 

= \ J {lnr 2 S T ab ^ }Sf dtAdr . (2.11) 

If the matter Lagrangian depends on the metric and its first derivatives only, 
this implies that (k = 8irG) 



ttiO r T Q T^o-b ^9ab 

E / C M = k — S 1 — 



/3=0 , f = N,S,/3. (2.12) 

Inserting the parametrization (|2.4J) for the metric now yields the desired 
result: 

El Cm = K^SN-^T^-T't), (2.13) 



V° s Cm = k^-T\, (2.14) 



E^C M = n—ST r t . (2.15) 
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These relations, together with eqs. ( |2.7| ) - (|2.9|), establish our assertion and 
enable us to write the Einstein equations in the form 

m' = -E°£ M , S' = E° m C M , mS = ^C M . (2.16) 

The fact that (3 may be set equal to zero after variation reflects the 
freedom to diagonalize the metric of the orbit space M. However, it is not 
surprising that (in a non-static situation) one looses information by using 
a diagonal metric in the effective action in the first place, that is, before 
performing the variation. 

Before we proceed, let us give two examples for which the above reasoning 
applies, since, in these cases, the matter action contains no derivatives of the 
metric. 



2.2 Examples 

As a first example we consider the Einstein- Yang-Mills (EYM) system. The 
matter action is 

S M = ^Tr J F A*F , (2.17) 

where e is the (dimensionless) gauge coupling, F = dA + A A A is the field 
strength assigned to the gauge potential A and a star * denotes the Hodge 
dual with respect to the spacetime metric g. For simplicity, we restrict our- 
selves to the gauge group SU(2). (The generalization to arbitrary gauge 
groups is straightforward applying the equations derived in A spheri- 

cally symmetric gauge potential A can then be represented as 

A = a r r + (Re(w) — 1) jr^, di? — Tasini? df\ + Im(w) {t# d$ + Tasini? d<^| , 

(2.18) 

where a is a one-form on M, w is & (complex) function on M and r r , r#, t v 
denote the spherical generators of SU(2), normalized such that [r r , r^] = t v . 
A gauge transformation with U = exp(xr r ), where \ is a function on M, 
preserves the form of the potential A and induces the transformations 

a^a + dx, w^e lx w. (2-19) 

Hence, a can be considered a gauge potential and w a Higgs field on M. 



With the parametrization (|2.18|) for the potential A one finds for the matter 
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action Sm, 



S 



Air 



M 



-(da\da) + (Dw\Dw) + 



(\w\ 



2r 2 



V 



(2.20) 



where we have introduced the covariant derivative Dw = dw — iaw. Using the 
parametrization ( |2.4| ) for the metric on M and adopting the temporal gauge 
a = a r dr, we now easily obtain the following expression for the effective 
matter Lagrangian C M : 



with 



M 



£ (0) + L 



M > 



(2.21) 



1 AO) 
or 



1 

S I N 



w\ + —a 



■2 



w 



ia r w 



2 (\w\ 2 

+ ■ 



2r 2 



2 (3NSRe{w(w' + ia r w)} 



(2.22) 
(2.23) 



and a 2 = AnG/e 2 . The Einstein equations are obtained from eq. (12.161). 
Applying the Euler-Lagrange operators E*, E ar on the diagonal part ?foi 
the matter Lagrangian also yields the YM equations: E/£j$ = E°£ M = for 
f = w,a r . (In addition to these equations one has the YM Gauss constraint, 
which got lost as a consequence of the gauge fixing.) 

As a second example we consider the Einstein-Skyrme model, for which 
the matter action Sm is (see and references therein) 



S 



M 



P 



Tr / A A *A 



1 



16e 2 



Tr / dA A *dA . 



(2.24) 



Here / and e are coupling constants, A = U 1 dU and U is a SU (2)-valued 
function on M (describing the pion field). The "hedgehog ansatz" 



U 



cos x — 2 r r sin x , 



(2.25) 



with a function x on the orbit space M, gives for the effective matter La- 
grangian C M , 

C-m = C-m + ' (2.26) 
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where 

= 4(^«X 2 ) -S(Nu X ' 2 + v), (2.27) 

- 2 C { $ = 2(3NSu X x', (2.28) 
or 

with a 2 = 2nG/e 2 and 

u = (/e)V + 2sin 2 x, u = (2(/e) V + sin 2 xf~^ • (2.29) 

The entire set of field equations (including the Skyrme equation) can be 
derived quickly from the total effective Lagrangian Co + Cm- 

In the examples above, the diagonal part C^ of the Lagrangian consists of 
a term proportional to S^ 1 and a term proportional to S, which, henceforth, 
will be called the "kinetic" and the "potential" parts, respectively. Note 
that the kinetic part is quadratic in the time derivatives of the matter fields 
and that it only contains terms proportional to S^ 1 or (iVS*) -1 , whereas the 
potential part consists only of terms proportional to S or (NS). 

2.3 Basic equations 

Motivated by the previous examples, we assume in the following that the 
diagonal part C := of the effective matter Lagrangian has the form 

C = ^{*(°> + ^*«} " S{NU + P) ' (2 - 30) 

with real matter fields X(o) an d X(i) an d positive definite quadratic forms x\o) 
and x\iy 

xfj) = ( X(j) , B (j)X(j) ) , 3 = 0, 1 , (2.31) 
where the inner products for the matter fields are denoted by ( • , • ). (For 
the EYM system we have, for instance, X(o) = a r, X(i) = (Rew,Imw) and 
5(o) = a 2 r 2 , _B(x) = 2a 2 diag(l, 1).) We further require (in view of the follow- 
ing sections) that the function P does not depend on spatial derivatives of 
the matter field x ■= (X(o),X(i)) : 

B {j) = B {j) ( X , X ',;r), U = U( X , x'\ r) , P = P( X ; r) . 
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With these assumptions for C it is possible to reconstruct the total effec- 
tive Lagrangian Cq + Cm up to first order in (3. We will show that 

C G + C M = S(m> - Pm) + C - (3( X , C x , } + 0{(3 2 ) , (2.32) 

where x '■= (X(o)j X(i)) an d C x > := dC/dx' ■ For the field equations E° {£ G + 
Cm} = 0, / = m, S, (3, x, this gives 

m' = -C s , (2.33) 

S' = C m , (2.34) 

m = -S^ix,/^), (2.35) 

(C X J = -(£ x ,y + £ x , (2.36) 

involving only the diagonal part C = C$ of the effective matter Lagrangian. 

To prove our claim, we show that the m equation ( |2.35| ) follows from 
the remaining ones. We do so by rewriting the system in Hamiltonian form. 
Defining the conjugate momenta tc := (7T(o), 7rm) by 

^(o) = C X{0) = £-B(o)X(o) , 7T ( i) = C X{1) = J^g B WX(i) , (2.37) 
the effective matter Hamiltonian is obtained by a Legendre transformation: 
W = (tt, x) ~ C , 

= tL^+N^+SiNU + P), (2.38) 
where tt? \, 7T?-n denote the quadratic forms 

^H 7 ^)'^) 7 ^))' J =0,1- (2.39) 

Observing that 

C m = -H m , C s = -U s = S-'H (2.40) 

and 

X = H 7T , C x = —H x , C x i = —H x > , (2-41) 
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) are equivalent to 


m = 




S' = 




m = 


S~ 1 (H n , H x > ) , 


X = 


'Htt , 


7T = 


(n x >y - n x . 



(2.42) 
(2.43) 
(2.44) 
(2.45) 
(2.46) 

Taking advantage of these formulae it is now rather simple to complete the 
proof. Using the Einstein equations ( |2.42| ), ( |2.43j ) and, subsequently, the 
matter equations ( |2.45|) , ( [2.46|) , we obtain 



(msy = (s^nys - mn ri 
= (x,n x ,y, 

which is the spatial derivative of eq. (|2.35|) . (It would have been considerably 
less convenient to use the Lagrangian formulation and it would have been 
even less so to work with the matter equations for an explicitly given model.) 

To summarize, the diagonal part ( |2.30|) of the matter Lagrangian uniquely 
determines the field equations ( |2.42j ) - ( |2.46[ ). The last four of these equations 
form a complete set of differential equations for m, S, ir, x, and, as is easily 
seen, imply that the constraint equation ( |2.42j ) (the Hamilton constraint) 
propagates. 

It will become clear in the following sections, that the above setup is 
also well adapted for treating linear perturbations of spherically symmetric 
selfgravitating configurations. 



3 The Pulsation Equation 

In this section we linearize the field equations (|2.33|) - ( |2.36|) in the vicinity of 
a static solution. As we will see, the only metric perturbations which enter 
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the first order matter equations are the variations of m and Ar, where A 
denotes the Laplace operator for the orbit space (M, g) . Using ( |2.33|) and 
fl2.34p one easily finds 



Ar=i(JV5)' = ^-i(4,+2P). (3.1) 



Together with the linearized Einstein equations, these metric perturbations 
can then be expressed in terms of the equilibrium solution and the pertur- 
bations of the matter fields. This enables one to derive a pulsation equation 
for the matter fields only. 

In the following, N, S, x, etc. refer to a static solution of ( |2.33| ) - ( |2.36| ), 
whereas time-dependent perturbations are denoted by SN, SS, 6x, etc. . 

Let us first discuss the metric perturbations 5m and 5Ar. Since both m 
and x are quantities of first order, the m equation ( |2.35| ) yields 

(S6mJ= -(C x/ ,5 X y. (3.2) 
On the other hand, the Einstein equations ( |2.33| ) and ( j2.34j ) imply that 

{S5m)' = C N 5N - S5C S 
= S (N5U + 5P) . 
Using the matter equation ( |2.36|) for the equilibrium solution, we see that 

(S5my = -(£ x ,,5xY. (3.3) 

Hence, 

d(S5m) = -d(£ x , , 5 X ) ■ (3.4) 

By assumption, P does not depend on derivatives of the matter fields. Equa- 
tion ( |3.4p thus yields 

6m = N(U x ,,6 X ) ■ (3.5) 

Here the integration constant has been set equal to zero. In the soliton case 
this is a consequence of the regularity requirement for the center. For black 
holes, this is the correct choice since we restrict our attention to variations 
with fixed position of the horizon. (The relation between eq. ( |3.4|) and the 
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first law for black holes is disussed in appendix B.) The variation of eq. ( |3 . 1| ) 
is now immediately found: 

2 2 
5Ar = — 5m — - SP 

2N 2 
= -pr(U*,8 X )--(P x ,6 X ). (3.6) 

Next, let us consider the matter equation fl2.36j ). Since, as will be shown 
below, the metric perturbations can be eliminated from the linearized matter 
equation, the latter can be written in the form 

T6 X + {U M + Ug}S X = 0, (3.7) 

where we have introduced the operators 

T6 X = NS5(£ X J, (3.8) 

Um5 X = NS6 X {(£ X ,)'-£ X }, (3.9) 

U G 6 X = NS6 g {(£ x ,)'-C x } . (3.10) 

Here, 5 X and 5 g denote the variations with respect to the matter fields x an d 
the metric g, respectively. The first two operators are immediately obtained 
from the definitions ( |3.8j ) and (|3.9|): 

T ^ ( N t m B° m ) ■ (UD 

Um = P*C/ X 'x'P* - t [P*, N S U x , x ] + NS 2 (NU XX + P xx ), (3.12) 

with the differential operator p,,, = —iNSd/dr. For the operator Ug we find 
from eq. ( |3.10| ) 

V G 5 X = NS 2 6 g {±{(£ x ,y-£ x }}, 

= -NS 2 {{U x ,)' - U x ) 5N - NS 2 U x > 5Ar . (3.13) 
The curly bracket in the last expression for Uc §X can be further simplified: 

N{(U X ,)' - U x } = (N + 2P-l)- U x , + P x , (3.14) 
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as can be seen from eqs. ( |2.36| ) and ( |3.1| ) for the unperturbed solution. 
Eventually we can use the expressions ( |3.5|) and (|3.6|) for 5m and 5Ar to 
obtain the operator XJ G from eq. ( |3.13| ), 



U 



2NS 2 



G 



[u x ,(P x , .) + P x (U x ,, .) + (2P-l)±Ux(U yf , •>}■ 



(3.15) 



Here we do not discuss mathematical properties, such as domains of def- 
inition or (essential) self-adjointness for the operators in eq. ( |3.7| ). This was 
done by some of us for the EYM system in a previous publication [|16 . 



4 Mass Variation 

In this section we demonstrate that the second variation of the total mass 
yields the same expression for the fluctuation operator as the one previously 
derived by means of linearizing the field equations. First of all, we recall that 
the Komar formula provides one with an expression for the total mass M in 
terms of the (asymptotically) timelike Killing field k. Defining the "local" 
mass M(r) by the Komar expression over a 2-sphere with coordinate radius 
r (0), 

M(r) = [*dk = ^(NS 2 )' = mS + r 2 NS f - rSrri , (4.1) 

one has M = linioo M(r) = lim^ mS. (Note that asymptotic flatness implies 
linioo r 2 S' = linioo rm' = 0). Thus, in a gauge where lim^ S(r) = 1, we 
obtain for configurations with 5x{oo) = 

5M = \imS5m = lim SN( U x > , 5 X ) = 0, (4.2) 

where we have also used eq. (|3.5| ) with 8m(0) = in the soliton, and 
8m{r H ) = in the black hole case. Since 5M = for static solutions, 
we obtain the second order formula 

1 1 r°° 

/•; M ■ -d-.\l M ■ - (S5 2 m)'dr, (4.3) 

2 2 Jro 

where r = for solitons and r = r# for black holes. 
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Our aim is to establish the relation 

q {S5 2 m)'dr = J {(6 X ,T5 X ) + (5 X , (U M + U G )5 X ) } , (4.4) 

where the operators T and XJm + Ug in the kinetic and the potential parts of 
the fluctuation operator are the same as in formula (|3.7|) . This is most easily 
achieved by applying the Hamiltonian formulation presented above. Using 
the Hamiltonian (|2.38| ), the m' equation (|2.42|) becomes 



m 



Unf 0)+ Nn 2 {1) } + NU + P. (4.5) 



Since x is of first order, we have Sn^ = S 1 B^5x(o) and Sn^ = (SN) 
Sx(i), from which we obtain 

S 2 m ' = S- 2 ( 5 X{ o) , B {0) 5x(f» ) + S- 2 ( 5x { i) , N~ l B {l) 5 X{ i) ) + S 2 (NU + P) . 

Recalling the definition Q3.11J) for T and using eq. (|2.34|) for S', S f = 2SU/r, 
and 5 2 m = —(r/2)5 2 N, this already yields the kinetic term in the expression 
( |4.4|) for the second variation of m, 

(S5 2 m)' = -^(S X ,T8 X ) + S{N5 2 U + 5 2 P) + 2S5N5U. (4.6) 

Note that the terms involving the second variation of N cancel, since 

(S5 2 m)' = -SU5 2 N + S{5 2 m)' . 

It remains to show that the integrals of the second and third term in this 
formula coincide with the corresponding expressions in eq. (|4.4|) . Using the 
definition (|3.9|) for Um, the second term can be rewritten as follows: 



S(N6 2 U + 5 2 P) = -5 2 X C 

= -(6 x £ x ,,5xY+ (5 X {(£ X ,)' - £ x } , 6 X ) 



(5 X ,U M 8 X ). 



NS 

Here and in the following " = " stands for equal up to terms whose spatial 
integration vanishes. (The integral over the last term in the first line of the 
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above formula does not contribute, since the variations of the matter fields 
are required to vanish at r = r$ and for r — > oo.) Hence, the integral of the 
second term in eq. (|4.6| ) assumes the desired form. 

Our final task is to establish the equivalence between the last terms in 
eqs. flOf ) and ( [4.4|) . We do so by writing SU = U x Sx + U x /5x' in order to 
integrate S5N5U by parts, 

SSNSU = (5 X ,{U X -(U X ,)})S5N - (6 X , U x >) (SSN)' . 

The first term on the r.h.s. of this equation is already of the desired form 
(cf. the first term in eq. (|3.13|) ). In order to obtain the correct expression 
for the last term in eq. ( }4.6| ), we have to add an additional S5U5N to the 
above formula. Hence, it remains to compute S5U5N — ( 5\ , U x > )(S5N)'. 
Using again S' = 2SU/r and eq. (|3.5|) to express SN in terms of Sx, we find 



S5U5N - (S X , U X/ )(S5N)' = -(S X , U x >) {nSS(2U /r) + (S5N)'} 

= -(8 X ,U X ,) [NS5{S'/S) + {SSN)'} 

= -(S X ,U X ,)S{(NS)'/S}S 

= -(S X ,U x ,)SArS , 

which is the desired expression, corresponding to the second term in eq. 



(|3.13|) . Thus, the preceding two formulae together imply 



2 SSUSN = S(S X , {(U x ,)' - U x } SN - NS 2 U x , SAr ) 

= j^(Sx,U G S X ), (4.7) 

where we have also used the expression ( ^.13| ) for the operator XJq. This 
completes the proof of the variation formula (|4.4j) . 

For the examples in section 2.2, with Lagrangians ( |2.22| ) and ( p.27[ ), one 
can immediately read off the operators T, Um and Ug, given in (|3.11|) , ( p,12j ) 
and (|3.15|) , and thus write down the explicit form of the pulsation equations 



as well as the second variation of the energy (^4). The results and 
their discussion can be found in |TB| and 
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Appendix 

A The gravitational Lagrangian 

This appendix is devoted to the dimensional reduction of the spherically 
symmetric Einstein-Hilbert action. Our aim is to derive the expression (|2.2| ) 
for the effective Lagrangian in terms of the intrinsically defined function 
r: M —>■ R. Recall that M = M/SO(3) is the orbit space with induced 
metric g, and spacetime (M, g) is a warped product of M and S 2 with metric 
g = g + r 2 g, where g = dVL 2 . In what follows, we will use a tilde for spacetime 
quantities and a hat for quantities on S 2 . 

A natural basis on the 2-dimensional orbit space is provided by the 1- 
forms dr and *dr, which are used to introduce the orthonormal diade {9°, 9 1 }, 

0° = -*6\ 9 1 = N' 1/2 dr, N = {dr\dr), (Al) 

where g = —9°®9° + 9 1 ®9 l . Here * and ( | ) denote the Hodge dual and the 
inner product with respect to g, respectively. Using Arrj = *Ar = d * dr, 
we immediately find 

d9° = - i- I -dN - Ar dr } A 9° , d9 1 = \- dN A 9 1 . (A2) 

Comparing this with the structure equations and taking advantage of dN A 
9° = (*dN) A 9 1 and dN A 9 1 = (*dN) A 9°, enables one to read off the 
expression for the connection form u = uj\ = u\ of M, 

u = — * (-dN - Ardr). (A3) 
N v 2 ) v ) 
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Using N = 1 — 2m/ r, this also yields 

r m a (dmldr) -i .... 
dr A uj = \ — - Ar - - — ^ ^ 77 , (A4) 

which will be used below. 

The general expression for the Ricci scalar of a product manifold with 
warping function r (with norm N = (dr\dr)) is 

R = R + ^(1 -N) - - Ar. (A5) 

Since the second structure equation for a 2-dimensional pseudo-Riemannian 
manifold reduces to Q\ = doj\ = ^Rr], the Ricci scalar of (M, g) is obtained 
from 

Rr] = 2du. (A6) 
Writing the volume form as fj = r] A r 2 dQ we thus have 



Rfj = 2 { r 2 du + 2 (— - r Ar) 77 1 A . (A7) 



Since we are interested in an effective Lagrangian involving no second deriva- 
tives of the metric fields, we integrate the first term by parts and subsequently 
use the expression (|A4|) for drAu. The Einstein-Hilbert action then becomes 

lQnS G = [ Rfj = 2 / d(r 2 ujA5Q) + 4 / ry A dfl , (A8) 

JAf JAf iM A 

where the terms involving the Laplacian of r have canceled each other. 

In order to obtain the correct effective Lagrangian after dimensional re- 
duction, it remains to subtract a boundary term B, involving the trace of 
the extrinsic curvature, from the above formula. In terms of the connection 
forms of spacetime, B is given by (see, e.g., fi~9H ) 

16tt73 = / Cj^ A *(0" A 9") , (A9) 

JdM 

where DM = DM x S 2 and it is assumed that the induced metric on dM is 
kept fixed. Using the connection forms 

= Co A B = co A B , Co A = r -f9 A , (A10) 



17 



and the expressions *(6 a A 6 h ) = -e ab r 2 dQ 2 , *(8 A A 6 B ) = e A BV and *(6 a A 
9 B ) = - * 9 a A *9 B , we obtain 

A *(0" A 0") = 2r 2 cu Arffi 2 - 4r * rfr A rffi 2 + 2^ 23 A 77 . (All) 

The last 3-form in the above expression does not contribute to the boundary 
integral, since it contains no volume-form dQ. In addition, the second term 
needs not be taken into account neither, since we do not consider variations 
with respect to r. (If one allows for variations with respect to r, this term 
contributes to the effective Lagrangian. As a matter of fact, in this way one 
obtains an effective Lagrangian which, in addition to the (afe)-components, 
yields also the angular components of the Einstein tensor.) 

Hence, considering coordinate, only the first term in eq. ( |A11| ) 

contributes under variations, 

l6n5B = 25[ r 2 uj A dVt . (A12) 

J dM 

Applying Stokes' theorem, we observe that this term cancels the first term 
in the variation of the dimensionally reduced Einstein-Hilbert action ( |A8|) . 
We thus obtain the desired result 

16^(5 - B) = 45 [ r, A dSl = I6n6 [ ^ dr \ , (A13) 

v ' Jm (dr\dr) Jm N 



completing the derivation of eq. (|2.2|) 



B The First Law 

In this appendix we establish that the first law for spherically symmetric, 
static black hole configurations reads 

1 r°° 

SM - —kSA = / (SSmYdr, (Bl) 

where k and A denote the surface gravity and the area of the horizon, respec- 
tively. This expression is valid for arbitrary matter models. For the theories 
under consideration in this article the term on the r.h.s. generically does not 
contribute. This is due to the fact that the integral is, by (|3.5|), equal to 
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lim 00 [S'A^( U x i , S\)], which usually vanishes as a consequence of asymptotic 
flatness. The above formula can, of course, be obtained from evaluating the 
general mass variation formula 



SM - -UcL4 — t~ f GTSgp, *k-^-S f *G(k) (B2) 

877 107T JT, 07T JT, 

(with G{k)^ = G^ u k v ) in the spherically symmetric metric used throughout 
this article. Here we shall, however, present a direct derivation of eq. ([B~T|), 
which is adapted to the spherical symmetry. 

Using the expression (11) for the "local" Komar mass and the require- 
ment of asymptotic flatness, M = linioo M(r) = lim 00 (m5'), we have 

POO 

6M = \im{SSm)= / (SSm)' dr + (S6m)(r H ) , (B3) 

where we have also used linioo S = 1 and linioo 8m = 5 linioo "m. Hence, in 
order to establish the desired result, we have to show that 

— k8A = (S5m){r H ). (B4) 

077 

To see this, we first note that 2m(rg) = yields 

(5m)(r H ) = i(l - 2m'{r H )) 5r H . (B5) 

In order to complete the derivation, we recall the general result M# = ^ft/l 
for the evaluation of the Komar integral over the Horizon. Comparing this 
with the expression Mjj = M(rjj) — [S(m — rm')](rjy) which is obtained 
from eq. Q4.1|) , the surface gravity of the horizon becomes 

'l-2m\r H )). (B6) 



2r 
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Using this and 5A = 8nr H 8r H in eq. ( p5|) yields the desired result 
Together with eq. (|B3|), this eventually establishes the variation formula (|B1|). 
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